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SUMMARY 


An  empirical  approach  is  used  lo  correlate  the  mixing  of 
isothermal  incompressible  coaxial  air  flows  in  straight  pipes  of  constant 
circular  cross-section.  An  assumed  decay  function  for  kinetic  momentum  is 
correlated  in  terms  of  initial  geometry  and  initial  velocities  of  the  mixing 
streams.  The  method  allows  the  rapid  estimation  of  quantities  of  engineering 
interest  including  the  axial  distribution  of  both  static  pressure  and  total 
momentum.  The  method  has  been  tested  against  published  data  on  thrust 
augmenting  ejectors  of  the  parietal  type,  and  found  to  provide  satisfactory 
reproduction  of  measured  performance  as  affected  by  mixing  pipe  length 
and  area  ratio. 

l 


RESUME 


On  exploite  line  approche  empirique  pour  etahlir  la  correlation  du 
melange  d’ocoulcmoiiis  d’air  coaxiaux  ineompressibles  et  isothermiques  dans 
des  cunduilcs  dmites  de  section  eireulaire  constante.  La  correlation,  etablic 
en  fonction  de  la  geometric  initiale  et  des  vi  tosses  initiales  lies  cotirants  de 
melange,  utilise  line  function  de  decroissaneo  hypothetique  de  la  quantile 
de  mouvement  cinetique.  La  method?  relenue  rend  possible  line  estimation 
rapide  des  grandeurs  d'mterel  technique.  v  compris  la  repartition axiak  de  la 
pressjon  statique  et  de  la  quantite  de  mouvement  totale.  Los  resultats 
obtonus  out  etc  compares  aux  donnees  publiees  sur  les  ejecteurs  d'augmen- 
tation  de  la  poussee  de  type  parietal,  et  on  a  trouve  que  cette  method? 
roprodiusait  de  lacon  sat isfaisante  la  performance  nesuree,  telle  que 


modifiee  par  la  longui  "r  de  melange  et  le  rapport  des  sections. 
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CONKINK.I)  MIMNC,  OK  COAXIAL  FLOWS 


1.0  INTRODUCTION 

Most  test  data  on  the  internal  aerodynamics  of  duct  components  relate  to  approach  flows 
with  thin  boundary  layers  or,  less  frequently,  fully  developed  velocity  profiles.  Little  information  is 
available  on  the  effects  of  entry  velocity  profiles  of  significant  non-uniformity  (referred  to  here  as 
'distorted'  velocity  profiles  and  characterized  by  values  of  momentum  coefficient.  >  1.02  say). 

Some  earlier  work  was  concerned  with  the  overall  effects  of  a  particular  type  of  entry  How 
distortion  on  the  static  pressure  recovery  of  diffusers  and  sudden  area  enlargements  (  Refs.  1  and  2 1. 
Subsequent  to  this,  it  appeared  useful  to  undertake  a  more  general  experimental  investigation  into  the 
behaviour  of  different  types  of  velocity  profile  distortion  in  the  flow  in  straight  pipes  of  circular  and 
annular  cross-section. 

The  process  of  confined  mixing  of  non-uniform  velocities  occurs  in  many  applications  of 
engineering  interest,  the  non-uniformities  arising  from  upstream  components  such  as  bends,  enlarge¬ 
ments.  diffusers,  manifolds,  turbomachinery,  ejectors,  etc.  Considerable  attention  has  been  directed 
in  the  past  toward  prediction  of  the  turbulent  mixing  process  by  relatively  sophisticated  calculation 
procedures.  I'hese  methods  necessarily  require  detailed  specification  of  the  flow  field  leaving  the 
distortion  generating  device  in  order  to  predict  the  parameters  of  interest  at  arbitrary  distances  down¬ 
stream.  It  is  noted,  however,  that  in  practice  appropriately  detailed  information  on  the  initial  flow 
field  is  not  generally  available  and  that  the  parameters  of  engineering  interest  in  the  mixing  duct  are 
frequently  confined  to  the  axial  distribution  of  static  pressure  and.  in  particular,  to  the  magnitude 
and  position  of  peak  static  pressure  recovery  together  with  some  indication  of  the  velocity  pvoiih 
characteristics  as  tins  region  is  approached. 

From  the'  standpoint  of  practical  need  therefore,  it  appeared  worthwhile  to  investigate  the 
possible  usefulness  of  an  empirical  correlation  technique  based  on  experimental  observations.  1  i i e 
method  assumes  a  plausible  decay  function  for  a  single  parameter  characterizing  the  velocity  profile 
and  Uses  the  momentum  equation  to  predict  the  longitudinal  static  pressure  distribution  in  the  mixing 
duet.  Appropriate  constants  are  derived  from  experimental  results  and  correlated  for  general  applica¬ 
bility. 


Attention  is  confined  here  to  mixing  of  coaxial  streams  in  a  dial  of  constant  circular  cross- 
section  under  essentially  incompressible  and  isothermal  conditions.  Sudden  enlargements  are  included 
as  a  special  case  where  one  of  the  coaxial  stream  velocities  is  zero.  <  filler  types  of  distorted  velocity 
profile,  and  the  effects  of  annular  rather  than  circular  mixing  pipes,  are  dealt  with  in  a  following 
report . 


2.0  ANALYSIS 

2.1  The  Momentum  Equation  in  Confined  Mixing 

The  momentum  equation  for  isothermal  incompressible  flow  between  planes  j  and  k  in  a 
straight  pipe  of  circular  cross-section  and  constant  area  may  lie  written: 


p>V-A  +  p,  A  =  p\pV  -  A  +  pkA  +  T-I)(Xk-  X,l 
where  ~  is  an  effective  value  of  wall  shear  stress  between  planes  j  and  k. 
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Th(*  mixing  of  coaxial  flows  of  differing  velocities  results  in  a  characteristic  type  of  longitu¬ 
dinal  static  pressure  distribution  which  is  illustrated  in  Figure  1  for  several  geometries  and  velocity 
distributions.  It  is  apparent  that  the  mixing  process  incorporates  three  main  regions,  illustrated  sche¬ 
matically  in  Figure  2.  An  initial  region,  characterized  by  relatively  minor  changes  in  static  pressure 
(and.  by  implication,  momentum  coefficient)  is  followed  by  a  region  of  rapid  mixing  in  which  the 
principal  rise  of  static  pressure  (and  corresponding  decay  of  kinetic  momentum)  occurs,  and  in  which 
a  maximum  value  of  (’  is  attained.  In  the  final  region  the  essentially  flat  and  highly  turbulent  velocity 
profile  near  till1  position  of  peak  static  pressure  recovery  readjusts  to  become  a  fully  developed  veloc¬ 
ity  profile  far  downstream.  In  the  vicinity  of  point  'a'  the  constant  gradient  of  static  pressure 
coefficient  is  associated  with  a  constant  friction  factor,  f.  and  a  momentum  coefficient,  p  .  corres¬ 
ponding  to  effectively  fully  developed  flow.  It  is  noted  that  p’  is  very  close  to  unity  Ip’  '  1.01  to  1.02 
depending  on  Reynolds  number).  If  friction  factor  were  constant  throughout  the  mixing  pipe  i at  f  I. 
the  datum  value  of  p  at  the  initial  mixing  plane  ( , J , ,  |  would  be  given  by  Relation  ( 1  I  as: 


-I dp  d’l  -  -dp  "  ,  +  lf  x., 


P,  •  (>! 

where  =  and  subscript  1  refers  to  '.he  initial  mixing  plane. 

1  2pY: 

In  practice  the  local  value  of  friction  factor  u  ill  differ  from  f  .  part  icularly  m  the  initial 
mixing  region.  The  actual  value  of  the  initial  momentum  coefficient.  fj (  .  can  tic  written 


•  i  • 

— p  i 


I  d  l 


where  represents  at,  ncrenicntal  pressure  coefficient  i  posit  i\ ,  . . r  i legal  i\ e  i  due  to  initial  wall  friction 
factor  differences  lr. >m  f  . 

If  it  is  a  -  -timed  that  these  differences  are  eon  fined  ,.d  n  glen .  i  e.  t  lie  w  all  friction 

factor  in  t  he  rapid  mixing  region  is  cssent  ially  constant  al  1  tie:  •  n.  ea,  mi  a  ntum  «  oeffi,  i  > i !  at  an 
arbitrary  pon  mi  m  'lie  rapid  mixing  region  may  i>i-  .  >1  a  a:  a  :  '  i  ■  e.  lb  latnni  i  1  i  ,e 


-(d>„  '  d  i  2p„  C(.  (  ,  •  if  ix  x.  . 
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where  ( ’ 

1  2pV: 


20„ 


201, 


-  If 


l  la  I 


or,  from  (2l: 
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I'll  is  expression  allows  the  evaluation  of  the  momentum  coefficient,  pn ,  at  any  point  in  the 
rapid  mixing  region  from  measurements  of  static  pressure  coefficient.  The  present  prediction  method 
is  based  on  fitting  the  experimentally  derived  values  of  tilt1  momentum  coefficient,  along  the  pipe, 
with  a  plausible  decay  function. 

2.2  The  Momentum  Coefficient  Decay  Model 

Assessment  of  the  measured  data  (e.g.  Fig.  1 1  suggests  that  momentum  decay  characteristics 
in  the  rapid  mixing  region  can  be  well  represented  by  an  expression  of  the  form: 


dP 

dx 


=  -kp 


Integrating  and  introducing  convenient  constants  gives 


p  —  p  |  e  ^  *  x  1 


(3) 


The  essence  of  the  prediction  method  lies  in  the  correlation  of  the  constants  k  and  d  with 
geometry  and  initial  velocity  profile.  In  turn  this  depends  on  the  adequacy  of  the  postulated  model 
and  the  degree  to  which  k  and  d  reflect  actual  physical  processes. 

In  the  initial  region  of  small  momentum  decay  it  is  assumed  for  consistency  and  complete¬ 
ness  that 

dP 

-  =  -zkp  where  z  «  1 

dx 

On  integrating  this  leads  to 

P  =  £Jje-/k'  (Hi 

d 

for  0  <  x  <  - 

1  —  z 


i.e.  d  is  closely  representative  of  the  initial  region  length. 

The  complete  model  is  summarized  in  Figure  3.  It  is  noted  that  if  k.  d  and  A  are  known 

then  C„  in  the  rapid  mixing  region,  for  given  (3,  ,  follows  from  Equations  (3),  (  la)  and  (5).  In  parti- 
1  n 

cular  the  position  and  magnitude  of  the  peak  static  pressure  is  given  by 
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x  =  d  +  —  tn 
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If  x  -  A 
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At  the  peak,  P  = 


2f 
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Tlu>  constants  d  and  k  were  evaluated  from  the  measured  data  by  fitting  each  faired  pressure 
rise  curve  at  two  points  in  the  rapid  mixing  region.  The  points  were  arbitrarily  selected  in  the  area  of 
greatest  interest,  at  5()b  a ml  90b  C  I  lien  k  and  d  followed  from  application  at  each  point  of 
F.quat  ions  (  1 )  and  (  f>|. 

2.3  Initial  Momentum  Coefficient,  p,  ,  for  Coaxial  Hows 

The  isothermal  incompressible  coaxial  flows  considered  here  are  described  by  the  diameter 
ratio.  I).  and  the  inlet  velocity  ratio  V,  ;Vtl. 

For  the  ideal  case  of  negligible  inner  pipe  wall  thickness  and  uniform  inlet  velocity  profiles 

(1  !)’  I  +  f)\V 

*P|  >i.i «-.ii  r 

[(  1  -  1):  i  +  I)" i  V, 

For  non-uniform  individual  entry  profiles 

p,,<  1  -  l): )  +  p,l):t  V,  \ 

( 1  iV  i  +  D- 1  \’,  V(,  I 

It  is  noted  that  for  turbulent  entry  flows  p  and  p’.  are  elose  to  unity,  (p1  -  1  .02  for  fully  developed 
turbulent  pipe  Howl  and  p,  is.  in  General,  not  significantly  different  from  (jj|  luU  (j. 

ip',  ,|  is  shown  plotted  against  I)  for  various  V,  V,,  in  Figure  I.  The  marked  red ue turn  a,- 
imiial  velocity  ratio  approaches  unity  is  evident,  and  is  further  illustrated  in  Figure  fi  which  presents 
the  maximum  attainable  value  of  ip,  i)lU.  l(  as  a  function  of  initial  velocity  ratio.  It  is  seen  that  values 
of  V  V,,  or  V„  \b  in  excess  of  about  0.1  will  generally  result  in  p,  values  too  low  for  significant 
pressure  recovery  through  downstream  mixing. 

3.0  EXPERIMENTAL  DETAILS 


, .  V  , )- 


3.1  Apparatus 

The  mixing  pipe  consisted  of  7-1  2  ft.  (2.3  mi  of  1.8  in.  (3.2  mm  i  wall  brass  tubing  of 
nominal  1.1).  3.75  in.  (0.5  uni  with  a  semi -circular  bcllmouth  of  3/8  in.  (9.5  mini  radius  at  entry. 

I  he  inner  surface  was  honed  smooth  resulting  in  a  final  1.1).  of  3.701  in.  (9.553  rail,  1  lie  pipe  was 
equipped  with  a  single  row  of  static  pressure  tappings  (Fig.  (ii.  and  provision  was  made  for  the 
mounting  of  traversing  gear  at  each  longitudinal  station.  Circumferential  traversing  and  measurement 
■/I'  ■'tat if  pressure  at  various  circumferential  positions  could  be  accomplished  by  pipe  rotation,  which 
vu.-  facilitated  by  roller  supports  and  a  clamped  slip  flange  with  an  ‘O'  ring  seal  at  the  downstream 
.unction.  Additional  static  pressure  tappings  were  provided  at  the  initial  mixing  plane,  located  21  I 
iianu  ter-  from  the  pipe  entrance.  The  pipe  was  sucked  from  atmosphere  through  a  b  in.  ( 15.2  ran 
■liameter  connection  containing  a  standard  venturi  meter  and  leading  via  a  gate  valve  for  control  to  a 
5  ft.  i  l(i. (>7  mi  diameter  continuously  evacuated  sphere.  Choked  flow  through  the  gate  valve  re- 
uuted  in  particularly  steady  flow  conditions. 

'  ariotis  coaxial  flow  situations  were  provided  by  a  series  of  inner  pipes  supplied  by  blown 
•  nr  a  id  terminating  at  the  initial  mixing  plane.  Outer  wall  surfaces  of  the  inner  pipes  were  machined 
town  to  give  wall  thicknesses  of  1  lb  in.  t  1  ,b  mm)  or  less  over  the  final  1  1  in.  I  35  uni  of  the  pipes, 
m  addition,  wakes  from  the  pipe  walls  were  minimized  by  tapering  the  outer  surface  tat  less  than  1  > 
.  i  achieve  a  final  wall  thickness  of  about  0.005  in.  (0.13  mm  I  at  the  end  of  each  pipe,  i  let  ails  of 
mier  mpe  dimensions  and  of  geometry  near  the  initial  mixing  plane  are  also  shown  m  Figure  b. 


I 


I 
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Mass  flows  in  the  inner  pipes  were  measured  by  a  standard  orifice  installation  well  upstream, 
and  a  subsequent  straight  run  of  inner  pipe  7  ft.  8  in.  (2.46  m,  29  to  90  dia.l  long  was  employed  in 
each  setup  to  ensure  essentially  fully  developed  flow  at  inner  pipe  exit.  Considerable  care  was  taken 
to  position  the  inner  and  mixing  pipes  in  a  truly  coaxial  fashion.  Boundary  layer  trip  strips  were 
located  on  the  outer  wall  of  tlu-  inner  pipe  and  on  the  mixing  pipe  bellmouth  as  shown  in  Figure  0. 

Air  temperatures  were  measured  by  thermocouples  adjacent  to  the  flow  measuring  device 
and  near  the  bellmouth  inlet.  With  the  exception  of  orifice  and  venturi  differentials,  which  were 
displayed  on  water  manometers,  pressures  were  measured  via  a  scanivalve,  a  0-1  psi  (0-7  kl'a)  differ¬ 
ential  transducer  and  digital  voltmeter. 

4.0  PROCEDURE  AND  RESULTS 
4.1  Results  with  V,  >  V„ 

Most  of  the  tests  were  conducted  at  a  constant  Reynolds  number  of  about  10'  in  the 
mixing  pipe.  i.e.  at  a  fixed  total  flow  rate.  Mean  velocity  in  the  mixing  pipe  was  usually  about 
(50  ft /sec.  (IS  m/s).  Inner  velocities  were  limited  to  essentially  incompressible  levels  (■  .400  ft  sec. 

(90  m,s)|.  The  suction  gate  valve  was  set  to  give  the  desired  overall  flow  rate.  The  velocity  ratio 
between  the  inner  and  outer  streams  could  then  be  varied  by  adjustment  of  the  air  supply  control 
valve  on  the  inner  pipe.  The  limiting  condition  of  V(1/V,  -  0,  corresponding  to  a  sudden  enlargement 
of  the  inner  pipe,  was  difficult  to  establish,  implying  as  it  did  equality  of  inner  and  mixing  pipe  How 
rates.  The  problem  was  circumvented  by  constructing  a  series  of  three  flanged  entry  pipes  I  details 
given  in  Fig.  (5)  to  investigate  the  sudden  enlargement  situation  directly. 

In  general,  for  each  inner  pipe  in  turn  the  air  flow  rate  was  varied  to  cover  the  useful  range 
of  initial  momentum  coefficient,  and  the  longitudinal  static  pressure  distribution  for  each  combina¬ 
tion  of  total  and  inner  pipe  flow  rates  was  measured.  Experimental  results,  plotted  in  the  form  C 
vs  x.  are  shown  for  each  pipe  arrangement  with  various  velocity  ratios  in  Figures  7(a)  to  if).  Lines 
have  been  faired  through  the  experimental  points  to  assist  in  interpolation.  Only  flow  conditions 
which  embody  a  central  jet  ( i.e.  V,  >  V.,  )  are  covered  in  these  figures.  Discussion  of  the  effects  of 
peripheral  jets  corresponding  to  V,,  >  V,  is  deferred  to  Section  1.4. 

It  is  evident  from  these  results  that  axial  static  pressure  gradients  tend  to  become  constant 
shortly  after  the  peak  pressure  region,  the  gradients  being  about  10D  higher  than  would  be  expected 
with  fully  developed  flow  in  a  smooth  pipe  at  the  same  Reynolds  number.  In  view  of  the  care  taken 
in  finishing  the  test  section  it  seems  likely  that  this  value  of  pressure  drop  per  unit  length  reflects  the 
higher  than  normal  turbulence  level  and  associated  flatter  than  fully  developed  velocity  profile  pre¬ 
vailing  after  the  main  turbulent  mixing  process.  Hot  wire  measurements  suggested  axial  turbulence 
levels  on  the  centre-line  of  '.V '■  to  ID.  and  pitot  tube  traverses  gave  a  momentum  coefficient  of  about 
1.010  fc.f.  1  .oD  and  1.020  for  fully  developed  flow  at  the  test  Reynolds  number.  Ref.  4)  at  the  final 
downstream  measuring  plane  some  20  dia.  from  the  initial  mixing  plane.  As  mentioned  earlier  some 
indication  of  the  variation  of  wall  friction  can  be  inferred  by  extrapolating  the  line  of  wall  friction 
pressure  drop  back  to  the  initial  mixing  plane,  and  comparing  the  intercept  ( 2/3, ,  )  with  the  value 
determined  from  velocity  ratio  and  geometry  (2^3,  i.  Values  of  p,  were  calculated  allowing  for  inner 
pipe  wall  thickness  effects  and  appropriate  momentum  coefficients  on  inner  and  outer  flows,  file 
value  of  jj'  was  taken  as  1 .0  1  as  suggested  by  traverses  20  dia.  downstream  of  the  initial  mixing  plane, 
i  point  ‘a‘.  Equation  (  2)  >.  Calculated  values  of  2(3.  are  shown  on  the  ordinates  of  Figures  T (  a  I  to  (fl 
•‘or  direct  comparison  with  the  intercept  values.  In  general  the  agreement  is  close,  indicating  that, 
or  practical  purposes.  A  may  be  taken  as  zero  (  Eq.  I  41 1  and  substantiating  the  hypothesis  of  effec- 
;ve|v  constant  wall  friction  in  the  mixing  region.  (See  later  discussion  under  1.4.1 

1.2  Correlation  of  Results  (V,  T  V()) 

\s  indicated  in  Section  2.2  the  assumed  exponential  decay  function  for  p  was  fitted  through 
oomis  at  50' ;  and  90D  C  on  faired  curves  through  the  data  of  Figures  7(ai  to  I  f).  providing  values 
■fd  and  k  for  each  test  condition. 
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Calculated  values  of  d  are  shown  plotted  against  velocity  ratio  for  various  diameter  ratios 
in  Figure  8.  (experimental  ranges  of  velocity  ratio  and  diameter  ratio  were  0  V../V,  ■'  0.12  and 

0.272  <  1)  <  0.819  respectively.  The  data  indicate  a  tendency  for  d  to  increase  with  VM/V(,  suggesting 
a  correlation  of  the  form 

d  =  ( function  of  geometry)  •(  inverse  power  of  ( 1  -  V t ,  / V , ) ) 

In  Reference  -1,  Abbott  and  Kline  present  data  on  flow  in  two-dimensional  channels  con¬ 
taining  a  backward  facing  step  and  conclude  that  the  reattachment  distance  for  flow  over  the  step 
appears  to  scale  with  step  height.  Although  the  physical  interpretation  of  d  in  the  present  work  is 
ill-defined  and  the  axisymmetric  geometry  would  be  expected  to  introduce  some  differences,  it 

»„  -  D, 

seemed  reasonable  to  attempt  a  correlation  in  terms  of -  .  A  logarithmic  plot  of  d(  1  -  V()  /V, ) 

— 

against  ( 1  -  D)  shown  in  Figure  9  suggested  the  following  correlation  for  d 

d  =  4.2(1  -  D)5/J(l  -  V0/Vi)_1  (10) 

This  relation  is  compared  with  the  data  in  Figure  8. 

It  seemed  probable  that  the  exponential  rate  coefficient,  k,  would  be  strongly  influenced 
by  the  velocity  difference  between  the  inner  and  outer  streams.  A  correlation  for  k  was  postulated 
in  terms  of  this  difference  normalized  by  the  higher  (inner)  velocity.  A  logarithmic  plot  of  k  in  these 
terms  is  presented  in  Figure  10.  It  appears  that  a  suitable  correlation  could  take  the  form 

k  =  (function  of  geometry)  •  (1  -  V0  /V t ) -  • 5 
where  the  function  of  geometry  is  relatively  weak.  The  expression  finally  adopted  was 

1  "I  l  lb 

k  =  0.75(1-  V0/Vj)2  5  - - —  (11) 

D-(l  -  D2 ) 

Expression  11  is  shown  relative  to  the  experimental  points  in  Figure  11. 

It  is  noted  that  for  many  applications  with  0.8  >  D  >  0.4  a  satisfactory  approximation  is 
given  by  the  simplified  expression 

k  =  (1-  V()/V,)2  5  (12) 

Although,  for  present  purposes,  detailed  static  pressure  behaviour  in  the  initial  region  is  not 
of  prime  significance,  an  estimate  is  included  for  completeness  using  Equation  (6)  with  z  =  0.05  for 
all  conditions. 

The  predicted  behaviour  of  static  pressure  coefficient  with  streamwise  distance,  based  on 
the  foregoing  correlations,  is  compared  to  the  measured  data  in  Figures  12(a)  to  (f)  for  all  test  confi¬ 
gurations  and  flow  conditions.  Agreement  in  the  rapid  mixing  region  and  in  particular  in  the  vicinity 
of  C‘p  is  good. 

A  quantitative  indication  of  the  capability  of  the  present  correlation  is  afforded  by 
comparing  the  predicted  values  of  Cp  and  x  with  those  observed  in  practice  (Figs.  13  and  14).  The 
results  confirm  that  the  suggested  correlation  procedure  reproduces  these  practically  important 
quantities  to  a  useful  degree  of  accuracy.  In  general  the  value  of  C(,  is  predicted  essentially  without 
error,  while  the  more  critical  x  is  predicted  within  half  a  pipe  diameter. 


As  a  chock  on  the  applicability  of  the  present  method,  measured  data  on  coaxial  mixing  by 
Ka/.insky  and  Brighton  (Ref.  5)  are  compared  with  the  predictions  of  the  present  correlation  in 
Figure  15.  It  is  noted  that  although  the  value  of  I)  of  0.11)7  is  well  outside  the  range  in  which  the 
present  experiments  were  conducted,  a  satisfactory  degree  of  agreement  is  evident. 

1.3  Results  with  V(1  >  V, 

For  this  series  of  tests  the  air  supply  to  the  centre  pipe  was  disconnected  and  flow  was 
induced  from  atmosphere,  the  value  of  Y,  being  controlled  by  a  valve  at  the  entry  far  upstream  of  the 
test  section.  In  other  respects  the  test  procedure  was  essentially  as  described  earlier.  Fxperimental 
ranges  of  velocity  ratio  and  diameter  ratio  were  0  -  V,/Y„  0.315  and  0.712  <  II  '  0.907  res  pec - 

t  ively . 

Fxperimental  results  for  conditions  involving  a  "peripheral  jet’  (V,,  .  ■  V  |  are  shown  in 
Figures  10(a)  to  (cl  with  lines  faired  through  the  data  points  to  assist  in  interpolation.  The  main 
features  of  the  mixing  process  appear  broadly  similar  to  those  evident  in  the  earlier  results,  the  extent 
of  the  initial  delay  region  increasing,  and  the  subsequent  rate  of  momentum  exchange  decreasing, 
witii  increasing  velocity  ratio,  V,  V  ,  (i.e.  with  inner  and  outer  stream  velocities  becoming  more 
similar).  It  is  evident,  however,  that  the  rate  of  mixing  is  in  general  more  rapid  with  a  peripheral  jot 
than  with  a  central  primary  flow  (i.e.  Y(  >  Y  I.  (This  difference  is  well  illustrated  by  data  presented 
in  Figure  17  which  were  obtained  from  total  pressure  traverse's  at  various  axial  stations  in  the  mixing 
pipe  for  peripheral  and  centra!  jot  situations  of  equal  area  (1)  =  0.707)  and  similar  initial  lower  to 
higher  velocity  rati  )  land  hence  similar  p't  ).  The  apparent  effect  of  wall  friction  in  assisting  velocity 
profile  flattening  in  one  case  while  delaying  it  in  the  other  is  marked  and  appears  to  persist  through¬ 
out  the  range  of  geometry  tested.) 

As  before,  extrapolation  of  the  downstream  line  of  constant  pressure  gradient  back  to  the 
initial  mixing  plane  allows  an  assessment  of  the  overall  incremental  friction  force  upstream  through 
comparison  of  the  intercepts  Vu  with  calculated  values  of  2p,  .  In  contrast  to  the  results  for  V,  -  V 
significant  force  increments  are  indicated.  As  would  be  expected  these  are  positive  (.'.  7-  Oi  and 
increase  with  increasing  I). 

It  is  noted  that  in  the  treatment  of  experimental  data  it  is  implicitly  assumed  that  wall  static 
pressure  measurements  are  entirely  representative  of  flow  section  mean  static  pressure.  Significant 
radial  static  pressure  gradients,  particularly  at  the  initial  mixing  plane,  could  account  for  the  indicated 
differences  between  ,j,  and  p’l,.  Traverses  of  a  specially  designed  miniature  static  pressure  probe, 
however,  failed  to  establish  the  existence  of  any  significant  static  pressure  variations  other  than 
apparent  effects  of  high  turbulence,  confined  to  the  shear  layer  between  the  inner  and  outer  flows. 
F.videnee  that  the  observed  differences.  A,  were  in  fact  due  mainly  to  excess  upstream  wall  friction 
force  j  over  that  which  would  have  prevailed  had  local  friction  coefficient  been  constant  at  the  down¬ 
stream  value  throughout  the  mixing  pipe)  was  provided  by  a  series  of  measurements  using  a  Preston 
tube  I  inner  diameter  0.0  10  in.  ( 1  mm)). 

A  selection  of  results  for  three  geometries  is  shown  in  Figure  18  in  terms  of  f/f  vs.  x.  These 
indicate  a  high  Initial  value  of  wall  friction  (  f/f'  at  the  initial  plane  '  ( V  v , .  VT-  )  with  a  subsequent 
rapid  drop  to  the  final  value.  The  assumption  that  increased  frictional  effects  are  largely  confined  to 
the  initial  region  appears  to  be  substantiated.  Estimates  of  the  overall  reduction  in  pressure  coefficient 
due  to  initial  excess  friction  were  derived  from  these  curves  by  integration.  These  are  shown  compared 
with  the  indicated  values  of  A,  Equation  (3),  in  Figure  19. 

Although  uncertain* ies  attach  to  the  derivation  of  both  sots  of  data,  the  agreement  is 
sufficiently  close  to  suggest  that  the  individual  discrepancies,  A.  derive  almost  wholly  from  initial 
wall  friction  effects. 


i 


All  values  of  A  ( including  those  derived  from  integrated  skin  frietionl  are  shown  plotted 
against  the  velocity  ratio  V„/V  in  figure  20.  A  satisfaetory  correlation  lor  the  geometries  tested  is 
afforded  hy 


A  =  <>.5<V,,/V  li  •  (If  I 


It  is  noted  that 


1  -  V,/V„ 

<V„/V  -  n  =  - 

i/d-  -  (i  -  v,;v„ i 


(ldi 


This  is  shown  plotted  against  Vj/V,,  for  the  j3ltlJN  conditions  (Fig.  5)  in  Figure  21. 

Assuming  the  A  correlation,  Expression  ( 1  3),  is  physically  based,  it  would  seem  reasonable 
to  expect  it  to  apply  also  to  cases  for  which  V,  >  V,,.  The  factor  I  V  V  -  1 )  is  shown  for  these  cases 
also  in  Figure  21 .  It  is  seen  that  the*  implied  values  of  A  are  negative  and  considerably  smaller  than  lor 
V  ,  >  V,,  particularly  at  high  values  of  p,  .  This  would  seem  to  point  up  the  difficulties  of  deriving 
these  values  experimentally  (by  intercept)  and  to  justify  their  (conservative)  neglect  for  practical 
purposes. 


The  overall  reduction  in  pressure  coefficient  corresponding  to  A  is  assumed  to  have  taken 
place  by  the  end  of  the  initial  region  (  Fig.  2).  i.e..  as  before,  wall  friction  coefficient  in  the  region  of 
prime  interest  (rapid  mixing.  Fig.  2)  is  assumed  constant  at  f  .  ( F.stimates  of  local  pressure  coefficient 
in  the  initial  region  require  some  sort  of  assumption  as  to  the  axial  distribution  of  this  reduction. 
Where,  for  completeness,  these  estimates  have  been  included  the  distribution  has  been  assumed  linear 
with  x,  over  0  <  x  <  t’.) 


1.1  Correlation  of  Results  (V,,  >  V, ) 

Curve  fitting  at  50b  and  90b  C(,  points,  as  before,  led  to  values  for  d  and  k  which  were 
correlated  by  a  procedure  similar  to  that  employed  with  the  results  for  V,  >  Y  A  logarithmic  plot 
of  d(  1  -  Vj/V,, )  against  the  effective  non-dimensional  step  height  I)  is  presented  in  Figure  22  and 
indicates  a  similar  relationship  to  that  observed  for  V,  >  V,,  (  Fig.  9).  Data  were  correlated,  at  least 
qualitatively,  by  the  expression 


d  =  1.1(D)5  c  (i  -  Vj/V„r' 


(14) 


which  is  shown  in  Figure  23  together  with  the  experimental  points. 

Experimentally  derived  values  of  k  are  shown  in  Figure  24  in  terms  of  ( 1  -  V,,  V  )  and 
again  exhibit  a  similar  relationship  to  that  found  earlier,  leading  to  a  correlation 


k  =  1.9(1  -  Vj/V,,  )--5 


1 

D-(l -  D- ) 


(15) 


The  factor  1.9  (c.f.  0.75  in  Expression  ( ID)  quantifies  the  previously  noted  increased  mixing  rate 
pertaining  to  a  ‘peripheral  jet’  (V,,  >  V,)  compared  with  a  ‘central  jet’  (V,  >  V„  )  of  the  same  lower 
to  higher  velocity  ratio  and  the  same  /3 1  (i.e.  the  same  total  to  ‘primary’  area  ratio). 
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Expression  (15)  is  shown  relative  to  the  experimental  data  in  Figure  25.  (As  before  it  is 
noted  that  the  indieated  dependence  of  k  on  geometry  is  weak  and  the  use  of  a  simplified  expression, 
viz. 

k  =  2.5(1  -  V,/V„)’> 

would,  in  most  applications,  result  in  no  significant  differences  in  prediction,  at  least  for  0.5  <  D  ■'  0.9). 

The  value  of  z  adopted  earlier  (viz.  0.05)  was  found  satisfactory  in  providing  an  approx¬ 
imate  representation  of  the  initial  region.  With  this  final  element  it  was  possible  to  produce  the 
predictions  of  static  pressure  distribution  shown  compared  with  experimental  data  in  Figures  20(a) 
to  (c).  Agreement  is  generally  good.  Predicted  and  measured  values  of  the  magnitude  and  axial 
location  of  peak  static  pressure  recovery,  based  on  Relations  (8)  and  (7),  are  compared  in  Figures  27 
and  28  respectively.  The  correlation  appears  satisfactory. 

5.0  SUMMARY  OF  PREDICTION  METHOD 

In  summary,  for  the  mixing  of  coaxial  air  flows  in  a  duct  of  constant  circular  cross-section 
under  essentially  isothermal  incompressible  conditions,  and  given  V|(  V(1,  D,,  D(1 : 

Calculate  (3,  from  Relation  (9)  or  (9a). 

Estimate  /3'  (  rx  1.01)  and  f'  from  Reynolds  number. 

Then  for  x  >  ft' 

Cp  =  2/3,  [l-e-k,*-d>]  -  4f'x  -  A 


/3  =  /3,  e~k<  x-cJ  > 


where  for  V,  >  V, 


r\ 


k  =  (1-  V„/V,)-  s  for  0.4  <  D  <  0.8 


U  =  4.2(1  -  V/V.r'U  D)5 /•' 


A  =  0  (for  practical  purposes,  but  see  Section  4.3) 


and  for  V„  >  V, 


k  =  1 


■9(1-  V,/V„^4_  1  _■ 

LD2(1  -  D-)  . 


k  =  2.5(1-  Vj/V,,)-  -5  for  0.5  <  D  <  0.9 


d  =  1.4(1-  Vj/V,,)-'  (D)5/.i 


A  =  6.5(  V()/V  -  1)  •  ( 4 f ' ) 

D-(l  -  Vj/V0  )  ] 

=  6.5  - - -  <1f) 

.1-02(1  -  Vj/V,,)  . 


(  Vi 

Note:  K  = - where  z  «  1.  In  the  present  report  predictions  for  the  initial  region  have  been 

1  -  z 

included  for  completeness  by  taking  z  =  0.05  and  using,  for  0  <  x  <  Q', 


Cp  =  2)3,  (1-  e-^M  -  4f'x  - 


P  -  P  |  e” ' k  N 


6.0  APPLICATION  OF  PREDICTION  METHOD  TO  EJECTOR  PERFORMANCE 

I  he  foregoing  correlations  for  mixing  of  coaxial  air  flows  may  be  applied,  for  example,  to 
the  estimation  of  axisymmctric  ejector  performance  under  incompressible  and  isothermal  conditions 
I  for  eases  with  parallel  mixing  pipes  and  either  central  or  parietal  primary  flow).  In  particular  they 
allow  estimates  of  the  effect  of  mixing  pipe  length.  By  way  of  illustration,  such  ejector  calculations 
have  been  carried  oul  and  compared,  where  possible,  with  experimental  results  from  another  source. 

for  simplicity  the  calculations  pertain  in  general  to  stationary  conditions  with  uniform 
entry  flows,  negligible  -n  oiidary  flow  entry  loss,  and  p'  =  1.0.  (Other  postulated  conditions  can  be  readily 
incorporated  i  For ..  given  mixing  pipe  length  L,  mixing  pipe  area/primary  area  ratio,  S,  and  down¬ 
stream  friction  la.  tor  !  .  the  calculations  proceed  in  principle  by  relating  entry  secondary  dynamic 
pressure  to  the  me  an  statu  pressure  at  the  initial  mixing  plane  through  the  identity  of  secondary 
entry  total  pressure  and  mixing  pipe  exit  static  pressure.  The  resulting  expression  for  initial  velocity 
ratio  contains  .  d  and  k  which  tire  themselves  related  to  velocity  ratio  through  the  correlations. 

An  iterative  procedure  is  therefore  required.  Starting  from  reasonable  initial  approximations. 
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t* .  A  =  0,  (1  =  0.  k  1,  convergence  is  generally  rapid  (.'t  or  1  iterations),  leading'  to  a  solution  for 
initial  veloeitv  ratio  anil  lienee  mass  flow  ratio.  The  final  values  of  d  and  k  allow  estimation  of  the 
momentum  eoeffieient  along  the  mixing  pipe,  in  particular  at  mixing  pipe  exit,  leading  to  thrust 
augmentation  ratio,  *1*  (ratio  of  ejector  thrust  to  the  thrust  corresponding  to  the  primary  flow  rate  at 
the  primary  total  pressure  discharging  directly  to  mixing  tube  exit  static  pressure). 

Some  sample  calculated  results  for  both  central  and  parietal  primary  jets  over  a  range  of 
values  of  mixing  pipe  length.  L.  at  a  typical  ejector  area  ratio.  S  ~  25.  are  shown  in  Figures  29  to  ,'il  . 
Figure  29  shows  the  calculated  values  of  k,  d.  A,  \’s;\  gj ,  .  p,  and  x  at  each  selected  value  of  1,. 

The  associated  performance  factors  <1>  and  'l'  are  shown  in  Figure  50.  It  is  noted  that  while  the  greater 
decay  rate  (k)  pertaining  to  the  parietal  arrangement  leads  to  a  considerably  shorter  mixing  pipe 
length  for  maximum  performance,  the  influence  of  increased  initial  friction  |  .\l  results  m  generally 
inferior  performance  particularly  with  respect  to  thrust.  The  differences  in  predicted  overall  perform¬ 
ance  between  the  two  ejector  type's  are  further  illustrated  in  Figure  51  where  ti,e  same  data  are 
plotted  as  <I>  vs  'F  and  shown  relative  to  ideal  performance. 

The  results  of  calculations  for  other  values  of  area  ratio.  S.  are  shown  in  Figure  52.  For 
clarity  only  the  maximum  values  of  the  performance  factors  <l>  and  'F  (together  with  the  associated 
values  of  mixing  pipe  length)  are  presented  in  this  Figure.  The  curves  exhibit  the  >ame  qualitative 
differences  between  the  ejector  types  as  observed  for  S  25.  with  the  additional  result  that  the  pre¬ 
dicted  maximum  <l>  attainable  with  the  parietal  type  declines  with  increase  in  area  ratio  beyond 
S  -  10. 


A  literature  survey  revealed  only  one  source  of  parietal  ejector  performance  under  appro¬ 
priately  isothermal,  incompressible  conditions  and  with  varying  mixing  pipe  length  I  Reid.  Ref.  <>). 
Experimental  data  from  this  reference  are  shown  compart'd  with  calculations  using  the  present 
method  in  Figures  34  and  35.  Figure  34  shows  thrust  augmentation  ratio.  'I*,  againt  I,  for  parietal 
ejectors  with  S  =  10  and  S  =  20.  The  generally  poor  level  of  thrust  performance  predicted  for  this 
type  of  ejector  is  well  substantiated  by  the  experimental  data.  (Calculated  values  of  <1*  are  shown 
both  for  uniform  secondary  entry  velocity.  12  =  1.0  and  for  a  more  realistic  secondary  entry  profile. 
S2  =  1.1.)  The  experimentally  observed  effect  of  mixing  pipe  length  is  closely  reproduced  by  the  pre¬ 
diction  method  as  is  the  adverse  effect  of  increasing  area  ratio  from  S  =  10  to  S  :  20. 

Measured  overall  performance.  <1*  vs.  *F.  for  a  wide  range  of  area  ratio.  S  -  10  to  131.  is 
shown  compared  with  the  predictions  of  the  present  method  in  Figure  35.  In  general  the  agreement 
is  surprisingly  good,  both  quantitatively  and  qualitatively,  particularly  in  view  of  the  extrapolations 
involved  in  applying  the  present  correlations  to  parietal  geometries  of  considerably  higher  values  of 
I)  (up  to  0.996)  than  those  used  in  the  present  experiments  (D  >  0.907). 


7.0  CONCLUSIONS 

An  empirical  approach  has  been  adopted  to  correlate  the  mixing  of  isothermal  incom¬ 
pressible  coaxial  air  flows  in  straight  pipes  of  constant  circular  cross-section.  An  assumed  decay 
function  for  kinetic  momentum  is  correlated  in  terms  of  initial  geometry  and  initial  velocities  of 
the  mixing  streams.  The  method  allows  the  rapid  estimation  of  quantities  of  engineering  interest 
including  the  axial  distribution  of  both  static  pressure  and  total  momentum.  The  method  has  been 
tested  against  published  data  on  thrust  augmenting  ejectors  of  the  parietal  type,  and  found  to  provide 
satisfactory  reproduction  of  measured  performance  as  affected  by  mixing  pipe  length  and  area  ratio. 
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FIG.  2:  MIXING  PROCESS  NOMENCLATURE 


FIG  4  IDEAL  MOMENTUM  COEFFICIENT  AS  A  FUNCTION  OF  GEOMETRY  AND 
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FIG.  7(c):  MEASURED  ST  AT 


0.273 


FIG.  7(e):  MEASURED  STATIC  PRESSURE  DISTRIBUTION  (Vt  >  VQ) 
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G.  11:  FINAL  CORRELATION  OF  DERIVED  k  VALUES  (Vj  >  VD) 
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FIG.  12(a):  PREDICTED  AND  MEASURED  VALUES  OF  Cp  COMPARED  {V-  >  VJ 


FIG.  12(b):  PREDICTED  AND  MEASURED  VALUES  OF  C  COMPARED  (V 
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FIG.  12(f):  PREDICTED  AND  MEASURED  VALUES  OF  C  COMPARED  (Vc  =  0) 
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FIG.  14:  PREDICTED  AND  MEASURED  VALUES  OF  x  COMPARED  (V,  >  V„) 


-  37  - 


FIG.  15:  DATA  OF  REFERENCE  5  COMPARED  WITH  PREDICTION  OF  PRESENT  METHOD 


FIG.  16(b):  MEASURED  STATIC  PRESSURE  DISTRIBUTION 
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FIG.  16(c).  MEASURED  STATIC  PRESSURE  DISTRIBUTION  (VQ 


.2  0.3  A  0.4  0.5 

(FROM  INTEGRATED  WALL  FRICTION  MEASUREMENTS) 

FIG.  19:  ALTERNATIVE  DERIVATIONS  OF  A  COMPARED 


FIG.  23:  FINAL  CORRELATION  OF  DERIVED  d  VALUES  (V0  >  V,) 


FIG.  25:  FINAL  CORRELATION  OF  DERIVED  k  VALUES  (V( 
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FIG.  26(a).  PREDICTED  AND  MEASURED  VALUES  OF  Cp  COMPARED  (Vc  >  V,) 


FIG.  26{b) :  PREDICTED  AND  MEASURED  VALUES  OF  Cp  COMPARED  (VG  >  V,) 


FIG.  26(c):  PREDICTED  AND  MEASURED  VALUES  OF  Cp  COMPARED  (V0  >  V*) 
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FIG.  27:  PREDICTED  AND  MEASURED  \ 
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FIG.  32:  PREDICTED  MAXIMUM  MASS  FLOW  RATIO  AND  REQUIRED 
LENGTH  vs  AREA  RATIO 


